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SUMNARY 


The interactions of spacecraft systems with the sur- 
rounding plasma environment have baen studied analytically 
for three cases of current interest. These studies include 
calculating the impact of spacecraft generated plasmas on 
the main power system of a baseline Solar Electric Propulsion 
Stage (SEPS ) , the modeling of the physics of the neutraliza- 
tion of an ion thruster beam by a plasma bridge, and examining 
the physical and electrical effects of orbital arabient plasmas 
on the operation of an electrostatically controlled membrane 
mirror. Xn order to perform these studies use was made of 
the NASA Charging Analyzer Program (NASCAP) , as well as 
several other computer models and analytical estimates. The 
main result of the SEPS study was to show how charge exchange 
ion expansion can create a conducting channel between the 
thrusters and the solar arrays. A new fluid-like model was 
able to predict plasma potentials and temperatures measured 
near the main beam of an ion thruster and in the vicinity of 
a hollow cathode neutralizer. Power losses due to plasma 
currents were shown to be substantial for several proposed 
electrostatic antenna designs. 
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1. INTRODUCailON 


This is the Final Report on Tasks 1, 2, and 3 of 
Contract NAS3-21762| '’Additional Application of the NASCAP 
Code". The work va$ performed by Systems, Science and Soft- 
ware between 26 March 1979 and 9 October 1980. 

Much of the material contained in this final report 
was published previously in topical, monthly, and conference 
reports. Produced separately under this contract are two 
other final reports. One describes NASCAP extensions in- 
cluding the new Low Earth O^rbit (LEO) model. The other pro- 
vides a detailed analysis of the SCATHA spacecraft. Other 
relevant NASCAP documents include "A Three Dimensional Dynamic 
Study of Electrostatic Charging in Materials," NASA CR-135256, 
and "NASCAP User's Manual," NASA CR-135259. 

The interaction of ion engine generated plasmas with 
spacecraft systems is a problem fundamental to the develop- 
ment of a Solar Electric P^ropulsion Stage (SEPS) . In Chapter 
2 we present a complex multidimensional analysis of charge 
exchange ion expansion and interaction with solar panels. 

The calculations indicate possible difficulties with the use 
of high voltages on the arrays. 

The major technical accomplishment under this contract 
has been the development of a workable theory of ion engine 
neutralization. In Chapter 3 and Appendix A, a fluid-like 
model of electron flow is described along with the results 
of two test calculations. The agreement with experiment is 
compelling evidence of fluid-like behavior. 

The last study, contained in Chapter 4, is detailed 
analytical investigation of environmental interactions of 
the electrostatic antenna baseline design. Included are 
spacecraft charging , electromechanical , parasitic current , 


and sputtering investigations. The primary area of importance 
again turns out to be power drain in low earth orbit due to 
current paths through the dense ambient plasma. 


2. PARASITIC CURRENT LOSSES DUE TO SOLAR ELECTRIC 
PROPULSION GENERATED PLASMAS 

2 . I INTRODUCTION 

Solar elactric propulsion is a leading candidate for 
lifting large structures from shuttle orbit to geostationary 
altitude. A previous study demonstrated that plasma pro- 
duced by charge-exchange reactions virithin the ion beam may 
dominate the ambient environment near the spacecraft. The 
ad hoc nature of the previous work made clear the need for a 
more accurate treatment of the expansion of the charge- 
exchange plasma and the resultant solar array power losses. 

The calculations presented here are an improvement over pre- 
vious work in that they contain a predictive model for the 
charge-exchange plasma expansion, and a fully three-dimensional 
model for the stirtcture of the plasma sheath around the solar 
array wing. 

For this study we have considered a 30 cm diameter 
mercury thruster producing 2 A of beam current and 25 mA of 
charge-exchange current. The charge-exchange ions are emit- 
ted radially, with energies of 5-10 eV, in the first thruster 

14 —3 

diameter of the beam. This gives a density of ^2.5 x lo m 

at the beam edge, which exceeds the LEO plasma density of 
10 12 —3 

10-10 m . Thus the charge-exchange plasma will dominate 
the ambient to a distance of several meters from the thruster 
in LEO, and over the entire spacecraft for substantially more 
tenuous environments. 

The conditions in this plasma are long collision 

3 —3 

lengths ('vlO m) and short Debye lengths ('\<10 m) . In this 

regime we assume that the electron plasma adjusts itself to 

maintain isothermal, quasi-neutral conditions. This implies 

a barometric law potential, 

<J) = eZnin/Tip) (2.1) 
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where 6 is the electron temperature (eV ) , n the local plasma 
density f and n^ the ambient plcsma density. This potential 
causes exp»>:ttBion and diffusion of the ion cloud as it emerges 
from the beam edge. Far downstream the electrostatic force 
become email, So that the ion density takes the form 
f(6)/r , where 6 is the polar angle relative to the beam 
and r the distance from the engine. 

The ion plasma is modeled hydrodynamically, i.e., 
we assume the ion density and velocity to be a well-defined 
function of position, and the ion thermal motion to be un- 
important. The motion then satisfies the equations of con- 
tinuity of mass and momentum: 

H--V.(pv) (2.2) 

a peE 

^ (pv^) . -VMv^py) + ^ (2.3) 

where subscript a denotes a Cartesian component, and E is 
the electric field. The challenge is to develop numerical 
methods capable of finding the steady-state solution to 
Eqs. (2. 1-2. 3) in the R-Z (cylindrical) geometry appropriate 
to the problem. 

2.2 NUMERICAL METHODS 

We solve the problem on the computational mesh shown 
in Figure 2.1. The physical space is located within the solid 
boundary; mesh points outside the physical space are compu- 
tationally convenient and serve to maintain boundary condi- 
tions. Velocities are defined at the points indicated by 
crosses, which form an evenly spaced r-z mesh: 

= r© + i Ar (2.4a) 

= Zq + j Az (2.4b) 
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Densities are defined at the centroids of the quadrilaterals 
formed by the crosses: 


.(P)2 . 1 r(v)2 «(v)2\ 

i I ^’^i * ^i+l I 


(P) . 1 /-(v) ^ (V) \ 

j 2 \*j * *j+l/ 


(2.5a) 


(2.5b) 


Mass flukes, f, are defined at the arrows in a fully *^upwind'' 


sense: 


''r (-r’ -*]'”) - k{\ * '-r 

arlv, (rl»>,»r>) -(rl:> - .1^') v, (r<-) >) 


{4”' - ^r’) ''z{4i 


i+l'*j 


(2.6b) 


. (r<-»/r<v)) p(r<^>,z<P;) v^(r ^ <P> ) 


i-1 if ^ 0 . 

i if < 0 


(2.6c) 


-pD = pH'” 


j-1 if V, > 0 


if < 0 
' z 


(2.6d) 


With above definitions, it is straightforward to timestep 
the integral form of Eq. (2.2), using as control volumes the 
rectangles formed by the four crosses (velocity definition 
points) Surrounding a density definition point. 






Equation (2.3) is then used to construct dv/dt. We 
use a finite difference method to evaluate the divergence 
term as follows: 

a. Define mass fluxes flowing between the crosses in 
Figure 2. If i.e., normal to the arrows. This is done by 
taking weighted averages of nearby already defined fluxes: 


i+1 j 


k»l Jl*j-1 


4ri 




i i+1 

E E 

k»i-l )l=j 




( 2 . 8 ) 


b. Define densities associated with velocity definition 
points: 


k-i-1 Jl=j-1 ' ' 

These densities are needed at both their previous and current 
(just updated) values. 



c. 


Define velocities at the arrows: 


d. 


where 


(V) _ 


i+1 


+ (^i+1 - ^i^^) y( 


,Zj 


(V) 


(V) ,(v) 


»Z 


j 


(V) . -Cv) 


■j+x ‘i 






Determine the change in yelouity by 


Q 


■ - ^®k("k’!k)Yk + '^ij ? 


ID 


^ij volume 



Sj^ are the surfaces bounding that volume, and is the force 
at Specifically 

«ij ' ^ - ■'i-i’) - *i-x) 

\ 

+ n(r^P>2 . 


V.^(t) » p| . (t-dt) V. . (t-dt)/p. . (t) +a..(t) dt 

'^Xj Xj '^Xj Xj '^XJ 


d t 


^(t) X®lt("k‘?n) -k ’ 


aij(t) 


. (e/m) 5(r!''>,af>) . 


These formulae are accurate to first order in dt and second 
order in (dr,dz) . 

e. Additionally, provision is made for some velocity 
smoothing at each timestep. 
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2 . 3 RESULTS 

As an example, we calculated the expansion of an ion 
plume with initial conditions similar to a case measured by 
Hughes Research Laboratory. The initial conditions were 

Emission energy (radial) = 10 eV. 

Initial radius r^^^ ■ .15 m. 

Current density e 

( 0 z < 0 

(.123 exp (-z/. 22) z > 0 

(total current ^ 25 mA) 

Electron temperature 6^ = 1 eV. 

Mass (Hg ion) » 3.34 x lO”^^ kg. 

12 -3 

Ambient density n^ » 10 m . 

It follows that 

^r^^^ ,z j = 3095 m/sec 

- 2.48 X 10^^ m"^ 

W = 5.5 V. 

The calculation was done in three phases: 

a. .15 m < r < 1. 15 m; Ar = 0.05 m 

-1.0 m < z < 1.0 m; Az = 0.10 m 

b. 1.0m < r < 5.0 m; Ar = 0.2 m 

-2.0 m < z < 2.0 m; Az = 0.2 m 

c. 2.0 m < r < 12 m; Ar = 0.5 m 
-4.0 m < z < 4 m; Az = 0.4 m 
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Initial conditions for the second and third phases were ob- 
tained by interpolating data from the previous phase and re- 
normalizing to retain the correct total current. At each 
phase^ the calculation was carried out until steady-state 
densities and velocities were reached. 

The results are shown in Figures 2.2 and 2.3. It is 
seen that the initially asymmetric expansion (Figure 2.2) be- 
comes roughly spherical by a radius of '^'l m from the beam. 
This is seen in Figure 2.3a, in which spherical (r + z ) 
scaling maintains constant arrow length, and in Figure 2.3b, 
which indicates radial contours of equal density. The den- 
sity (m”^) beyond 1 m from the engine is reasonably approxi- 
mated by 

p » 10^^ (sine)^^/(r^ z^) 

where 6 is the angle from the beam direction. A closer in- 
spection of Figure 2.3b indicates that the plume extends 
further upstream than downstream. This is attributable to 
"pressure blowoff" at the upstream plume edge near the engine 
where the density gradient is high;, 

The plasma parameters determined from the above model 
are used as input to a fully three-dimensional computer pro- 
gram designed to predict current collection by high voltages 
in low temperature, short Debye length plasmas. ^ ^ This 

model uses an analytic, nonlinear space charge formulation, 
correct in both Debye screening and thin sheath limits, to 
determine the electrostatic potential and the boundary of 
the plasma sheath. The model allows the plasma temperature 
and density to vary in space. Figure 2.4 gives a sample of 
electrostatic potential contours near a solar array Wing, 
illustrating the asymmetry caused by the charge-exchange 
plasma being predominantly on one side of the wing. By 



1 1 1 V 








Figure 2.4 


, Potential contours about a 2 kV solar array in the 
presence of the charge exchange plume. Note that 
the potential is more shielded above the wing (where 
the charge exchange plasma is located) than below. 




tracking electrons inward from the sheath boundary, good 
parasitic current estimates are obtained. Iterating on the 
two stages of the calculation allows non-local effects to be 
included. Additionally, a high resolution capability is 
available to compute the current distribution over a complex 
pattern of solar cells. 

Table 2.1 gives sample results for an 8 m k 30 m, 25 kw 

solar array for several orientations and configurations. In 

these calculations, the spacecraft body was held at plasma 

ground and the solar cells were positively biased. Parasitic 

currents to both sides of each 10 m section of the wing were 

separately monitored. Several trbnds are apparent: (1) The 

inboard section of the wing draws most of the current, even 

-2 

though it is at a low voltage. This is because of the r 
expansion of the charge-exchange plasma. (2) The current to 
the outboard section is similar to the center section current, 
due to the large end effect in the tenuous plasma. (3) When 
the beam is in the plane of the panel, an increased loss is 
caused by the array's intersecting che charge-exchange pan- 
cake. It is apparent from these results that with kilovolt 
biases power losses of 'V'lO percent or more are likely, even 
with only one engine in operation, and that ameliorative 
measures should focus on the inboard portion of the solar 
arrays. 
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3. ION ENGINE NEUTRALIZER STUDY 
3.1 PHYSICAL MODEL 

The purpose of the present investigation is to deter- 
mine if simple models for electron transport give qualita- 
tively correct predictions of the plasma properties in the 
vicinity of an ion thruster. 

within a meter or so of the ion thruster, the electron 
densities are in the range 

10® < n < 10^® cm"® 

and their velocity distribution is characterized by tempera- 
ture 0 between about one and ten electron volts. The Debye 
length 

Xjj a 700 cm 

is typically small compared to distances L over which there 
is a substantial variation of macroscopic plasma properties 
such as density; potential, and temperatures. On the other 
hand, the mean free path for pair collisions 

a 10^^ E®-/® 0^/®/n cm, E < 0 

for t^lectrons of energy E (eV) is typically long compared to 
L, so that the behavior of the plasma is controlled by col- 
lective rather than collisional effects. Since « L, the 
plasma is quasineutral, departures from neutrality amounting 
roughly to 

/ X 

&n/n 'x. y ^ 10*^ , 

the space around the vehicle is strongly shielded from surface 




potantials. This ia in contrast to the aituation that prevails 
in charging of spacecraft in geosynchronous orbit where effects 
of space charge are entirely negligible and potentials are 
determined as solutions of Laplace’s equation. 

Although collisionless f thruster-generated plasmas 
exhibit macroscopic behavior similar in many respects to that 
of a collisional plasma. Such behavior is perhaps not totally 
unexpected in view of the fact that in both non-equilibrium 
and equilibrium plasmas electrons are scattered by fluctuating 
electric fields. A primary difference between the equil- 
ibrium and non-equilibrium cases is in the magnitude of the 
fluctuating fields. 

Several investigators have measured properties of 
thruster generated plasmas. in the experiments of 
Ogaws/ « on cesium ion beams neutralized by electrons 

from a hot wire, measurements were made of the density, 
potential, and electron temperature in the beam plasma. The 
potential difference between the neutralizer wire and the 
plasma could be varied by changing the position of the wire, 
the large potential differences (electron injection voltages) 
occurring when the wire was completely withdrawn from the 
beam plasma. An important result of the Ogawa experiments 
was that over a wide range of conditions electron density 
n and plasma potential (p were well correlated by the bsro- 
metric law 

n(r) « const exp (q4> (r) /kT) (3.1) 

where T is the electron temperature. The approximate validity 
of the barometric law was further verified by Kaufman. 


20 


since the barontetric law ia a thermal equilibrium 
concept# it can be completely valid only if the plaama ia 
iaothermal * The plaama ia only approximately iaothermal# 
noticeable deviation occurring aa one proceeda from the beam 
axis beyond the beam edge into the plaama formed by ambient 
and charge exchange iona. Kaufman obaervea an electron temp- 
erature in the charge exchange plaama only about half that 
in the beam, Ogawa^^^ and Sellen^^^ obtained meaaurable 
temperature variationa in the beam plaama over aeveral tena 
of centimetera in the downatream direction from the accelerator 
grid. The largeat deviationa from the barometric law were 
obaerved for large injection potentiala (>10 volta) . Probe 
traces in such cases also indicated departures of the electron 
spectrum from a Maxwellian shape. 

Probe measurements of the plasma potential in the 
thruster beam were made in SERT II flight and ground test 
experiments. The measurements show strong variation of plasma 
potential across the beam boundary about 20 cm downstream 
from the thruster grids. Such results are difficult to 
explain on the basis of a barometric law relationship unless 
the electron temperature or density vatiation from beaun center 
to beam edge is much hig ^r than might be expected from other 
measurements made in similar configurations. We anticipate# 
however# that such is not the case and# instead# that the 
observed behavior should be explained in terms of the 
anomolous resistivity of the thruster generated plasma to 
the flow of electron current. Thus# the primary Objective of 
the following sections of this chapter is to determine the 
capability of sii^i<ple transport models to explain# at least 
qualitatively# the experimental results. 


We first give the precise kinetic equation for the 
electron distribution and the first few moment equations 
expressing conservation of change# momentum# and energy. 

From there# we state the approximations leading to the trans- 
port equation which we eventually solve. The method of 
solution and the results of calculations will be the subject 
of a later report. 

Exact Equ a tion 

In principle# a kinetic approach based on the Vlasov- 
Boltzmann equation fully describes the spacecraft generated 
plasma. The complexity of such an approach# however# makes 
it impractical as a basis for conducting multidimensional 
calculations of plasma behavior. Besides# except near 
sources and collecting surfaces# where the distribution 
functic>n may change markedly# one should be able to ade- 
quately describe the plasma in terms of certain average 
properties of the distribution# such as temperature# density, 
and particle and heat fluxes. 

Quite generally the state of the plasma can be speci- 
fied by the distribution function ?_(r#v#t) that characterize 
each particle component a, where f^(r,v,t)dr dv represents 
the number of. particles of species a in the six dimensional 

4» *4 4> 4 

volume element dr dv about the position r#v in phase space. 
The kinetic equations which describe the distribution are 








( 3 . 2 ) 


For particles of mass m_ and charge q_ in an electric field 

4 3 - 3 

E and a magnetic field fi the "smoothed” force on a particle 
is 
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(« * ¥) 


(3.3) 


The effects of collisions between particles is taken account 
of by the collision term denoted here by C_. Here we attempt 
to describe the plasma in terms of its density 


n ■ ff{v,v,t) dv , 


(3.4) 


mean velocity 

^ - i/f(r,y,t)v dv (3.5) 

and certain higher velocity moments. For convenience, 
we have omitted the particle species subscript a. The first 
three moments of the kinetic equation yield conservation 
equations for particles, momentum and energy, as summarized 
below: 


Conservation of Particles 
1^ + V • n^ = 0 
Conservation of Momentum 

1^ (ran^) + V • mn(?^ + V»f - qn ^ 

where 

S mv G d^ , 
f = J*mf (v^v") d^ * * 


(3.6) 


(3.7) 


(3.8) 

(3.9) 
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is the pressure tensor. 


p * nm 


( 3.1 


is the scalar pressure, 


i . nm (55 - ?) 


( 3.1 


is the stress tensor, and v** » v - ^ represents the deviation 
of particle velocity from the mean velocity. 


Conservation of Energy 


3 t 


V • P * qnE*V + + Q 


( 3.1 


where 


^nm + pj ^ + 1 *^ + q 


( 3.1 


is the total energy flux. 


q = nm 


( 3 . 1 ^ 


is the heat flux, and 

Q=/2!^Cd5 (3.1! 

So far, the equations are quite general and involve no 
assumption that the gas is collision dominated or retains a 
Maxwellian spectrum of velocities. Separate conservation 
equations may be written not only for different particle species 
but also for different groups or particles of the same charge 
and mass. Primary electrons, for example, with significant 
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streaming energies could be treated as distinct from the main 
electron population which is taken to have a Maxwellian distri- 
bution of velocities. For the present, however, and until 
experimental or theoretical considerations dictate otherwise, 
we shall consider electrons as a whole and that their distri- 
bution varies slowly in space. 

Approximation for Electron Gas 

Consider that the plasma is in a steady state and that 
quasi-neutrality pertains throughout the bulk plasma (that is, 
away from electrodes and collecting surfaces) . The electrons 
and ions each satisfy the particle continuity equation 

V • n^^^ * 0 (i * +,-) (3.16) 

with = n__ * n. The momentum equation simplifies consider- 
ably if the electron drift velocity ^ is small compared to 

2 h 

the random velocity <v' >- and if the velocity distribution 
is nearly isotropic. Then, in the absence of magnetic fields, 

Vp + enll = $ 

where R represents the collisional drag between ions and 
electrons. In a classical plasma dominated by collisions, 

$ is composed of a part proportional to the relative motion 
u = between electrons and ions, leading to plasma 

resistivity, and to a thermal part proportional to the gra- 
dient of electron temperature, which is frequently neglected. 

In this approximation, Eq. (3.7) becomes 

Vp + enl = n ne J (3.17) 
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where 5 net current density and the plasma resistivity 

n is related to the electron-ion collision frequency by 



-E -i- 

47T v^. 


( 3 . 


and 0 ) is the electron plasma frequency. If the plasma is 
non-resistive and isothermal, Eq. (3.17) yields the baromet- 
ric law, Eq. (3.1). In this sense, Eq. (3.17), or more 
generally the complete electroh momentum equation, may be 
regarded as the generalization of the barometric law. 


If the plasma is not collision dominated, randomization 
of electron velocities may still occur through the enhanced 
levels of fluctuating fields in the plasma, such as occur for 
electron two-stream instabilities, or electronsion instabil- 
ities of the ion-acoustic or Eunemann type. These mecha- 
nisms are probably effective in coupling neutralizer electrons 
into the bulk plasma and in equalizing electron and ion mean 
drift velocities. They are often approximated by introducing 
an effective collision frequency; e.g., in the case of the 
Bunemann mechanism, which is effective where the relative 
drift velocity between electrons and ions exceeds the mean 
thermal velocity of electrons 


'■eff 


for scattering of electrons with velocity v from Langmuir 
oscillations with phase velocity Vp 


'eff 
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where is the energy density in the wave. These expressions 
should be compared with the classical collision frequency 
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where is the Coulomb logarithm, v^, the electron thermal 
velocity and 




<<0 


is the fluctuation energy per particle in an equilibrium 
electron gas. 


The determination of electron temperatures in the 
plasma requires consideration of the energy balance equation, 
Eq. (3.12). Making the same approximation in the equation 
expressing conservation of energy that were made in the 
momentum equation, yield 


V.f' * qnE*^ 


with 


(3.19) 


F » I p^ + q (3.20) 

Here 5/2 pV is the enthalpy flux of the drifting electrons, 
q the macroscopic heat flux, and is related to the 
effective joule heating associated with the relative motion 
of electrons and ions. The quantity R appears also in the 
electron momentum equation; for a pla?-ima controlled by 
collective effects it should be approximated in the energy 
equation in the same manner as in the momentum equation. The 
heat flux, q, contains new featuresi Classically, q, contains 
two terms; one proportional to the relative drift velocity 
between electrons and ions, and the other proportional to 
the gradient of electron temperature. 
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For th<e initial calculations, the results of vi^hich 
will be reported later# we ignore the drift contributions 
to the energy flux# the electron-ion heating and assume 

that the heat flux is proportional to the temperature gradi- 
ent* The energy balance equation thus assumes the simple 
form 

_ _ 2 

V*KV0 + “**"^eff ^^e “ ^i^ * ® 

3.2 ION ENGINE NEUTRALIZER CODE 

The basic physics of an ion engine neutralizer model 
was presented in the preceding section. This physical model 
has been incorporated into a two-dimensional (R-Z) computer 
code# which is described below. Two sample calculations 
have been performed with this code: neutralization in a 

thruster similar to SERT II is discussed in Section 3.3.# 
and a high resolution calculation of the temperature and 
density profiles near a hollow cathode neutralizer is pre- 
sented in Section 3.4. 

The ion thruster model has been incorporated into a 
two-dimensional (R-Z) computer code following the block dia- 
gram shown in Figure 3.1. For this initial version the ion 
currents and densities were assumed known. In a later ver- 
sion# it would be possible to allow a multicomponent ion 
composition to be determined self-consistently with the tem- 
perature and potential. The code operates entirely in MKS 
units. 

The code has been run interactively, with all relevant 
information on disk file. As long as previous information 
exists on disk, the program may be entered from the two noted 
entry points as well as the beginning. For developmental 
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purposes# it was found convenient to "hard-wire" many features 
of a particular problem into the code# while others are 
prompted for input. A flexible on-line graphics program# 
which plots information on the disk file# has also been de- 
veloped. 

At present# the code assumes ion velocities everywhere 
and ion currents at the input boundary to be known. The code 
then calculates plasma densities such that y* (ny) « 0 is 
numerically satisfied. Typically, ion velocities are taken 
to be either purely axial or to be radial from a point source 
on the axis exterior to the mesh. As the code requires non- 
zero plasma density everywhere# a background density of "slow" 
ions may be added. It should be possible to handle multiple 
ion species with interconversion fairly easily. 

The neutralizer is assumed to be a ring at specified 
distance from the axis# emitting a current of electrons equal 
to the ion beam current. (In practice# this gives much too 
low a current density near the neutralizer. For this reason# 
we constructed, a modified, high-resolution version of the 
code, focusing on the neutralizer to determine the temperature 
and potential profile near the neutralizer orifice. ) The net 
current in the plasma is given by 

'?net * ^~ion ” -electron^ (3.21) 

= o (-y (|) + I Vp) , . C3.22) 

where o is the plasma conductivity and p = nkT is the electron 
pressure. [For and ekT-»-6 (constant) we find ({> - 0Jln n.] 

The code determines electrostatic potentials by solving 
~*~net ~ (See Appendix B. ) It is necessary to iterate 

between this equation and the temperature equation, since 
the pressure is a function of temperature. 
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The plasma temperature satisfies the equation 


,2 

7.(-kV0) * 

where k is the thermal conductivity and the riqht-hand Side 
represents the ohmic generation of heat. For this prelim- 
inary version, convective heat transport has been neglected. 
On the various boundary regions, either isothermal or insu- 
lating boundary conditions may be specified. In practice, 
since k * nK^e”"^, the equation actually solved is 

n t2 

-v*K^v(e") ~ 

(J 

For convenience, the transport coefficients a and 
are calculated by a single isolated subroutine. The conduc- 
tivity a may depend on both density and temperature, and k' 
on density only. The present version assumes a relaxation 
rate proportional to the plasma frequency: 

rr . n V2 e^ 1 

” m 8.98a 

where the parameter a is taken to be 0.51. By the classical 
Weideman-Franz law, 

3 k ^ 

< * f a (f) T 

3 

If we measure temperature in eV, k = e, so that k' * y a. 


(3.23) 


(3.24) 


(3.25) 


(3.26) 
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3.3 BEAM NEUTRALIZATION CALCULATION 

A calculation was performed for neutralisation of a 
0.25Af purely axial beam of 3 keV mercury ions in a constant 
density background plasma. The beam had a radius of 7 cm, 
and a hot-wire ring neutralizer was placed at a 17-cm radius. 
These parameters were chosen to model the SERT II ion thruster. 

The plasma density (Figure 3.2) had a peak (central) density 

14 -3 14 -3 

of 16.1 X 10 in , and an ancient density of 1.0 x 10 m . 

The electrostatic potentials are shown in Figures 3.3 and 
3.4. The strong fields at the beam edge (Figure 3.3) are 
dominated by the barometric law, while the axial field at 
the beam entrance is needed to satisfy the boundary oonditipn 
of stationary electrons and moving ions. Figure 3.4 (poten- 
tials at the thruster plane) shows the potential dip at the 
neutralizer position. The current-vector plots (Figures 3.5 
and 3.6) indicate that the beam is neutralized by electrons 
entering from the side. The beam is 50 percent neutralized 
at ~ 15 cm downstream from the thruster. The temperature 
profile (Figure 3.7) was calculated with insulating boundary 
conditions at the thruster. The ma,ximum temperature occurs 
at the beam entrance, where the heat generation is greatest. 

The peak temperature was 5 eV, compared with a 1 eV assumed 
background. 

3 . 4 HIGH RESOLUTION CALCULATION OF TEMPERATURE AND 
DENSITY NEAR A HOLLOW CATHODE NEUTRALIZER 

A version of the codes was constructed with .002 m 
radial resolution and .004 m axial resolution to treat the 
region in front of the "keeper” of a hollow cathode neutral- 
izer. The parameters were chosen to approximate the plume 

f X41 

mode experimental conditions of Siegfried and Wilbur'^ ^ 
(Figures 3.8 and 3.9). A current of 2 ampere was assumed 
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Figure 3.5. Current vectors for SERT II 

ion thruster model. The length 
of each arrow is weakly dep>endent 
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to flow through a 1 cm radius keeper. (Numerical difficul- 
ties, since resolved, prevented simulation of a 0.5 cm radius 
keeper.) Isothermal (1 eV) temperature boundary conditions 
were imposed at the keeper plane. The ions were assumed to 
flow radially from a source 1 cm behind the keeper, with a 
density of 10 m at the keeper (Figure 3.10). Our re- 
sults (Figures 3.11 and 3.12) predict, in agreement with ex- 
periment, a sharp rise in both temperature and potential, 
followed by a gradual drop as one proceeds downstream from 
the keeper. The numerical agreement is quite good, and would 
probably improve considerably if the calculations were done 
with a smaller keeper. Substituting insulating thermal 
boundary conditions at the keeper plane (Figures 3.13 and 
3.14) predicts far too much plasma heating. 
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4. ELECTROSTATIC ANTENNA MODEL 


4,1 INTRODUCTION 


The electrostatic antenna# or electrostatically con- 
trolled membrane mirror (ECCM) ^ ia one of several proposed 
systems (e.g, , solar power satellites) involving the use of 
extremely high voltages in space. Laboratory and space ex- 
periments, as well as theoretical calculations, have indi- 
cated that Interactions of such high voltages with the space 
environment can have severe adverse effects. For this reason, 
the Plasma Physics Group at Systems, Science and Software 
was asked to perform a study of this particular system with 
emphasis on identifying potontlal problems associated with 
spacecraft-environment interactions . 

At present, specifications for an ECCM spacecraft are 
in the concept state. Its construction, orbit, and even its 
use are largely undetermined. The first systems study is 
not yet under way. The results of our investigation indicate 
that the ECCM is a viable concept in that we have found no 
potential problem that cannot be circumvented (or at least 
substantially alleviated) by careful design. However, sput- 
tering in LEO is of uncertain magnitude and warrants further 
consideration. We feel that the information presented herein 
should be important input to any future systems or design 
studies. 




In Section 4.2, we summarize the basic concept of the 
ECCM, with particular emphasis on how various design features 
relate to spacecraft-environment interactions. Section 4.3 
describes the plasma environment in low-earth and geosyn- 
chronous orbit. Section 4.4 discusses spacecraft charging: 
the effects it may have on the membrane shape, as well as 
the possibility of actual damage. Section 4.5 discusses 
parasitic power losses involved in maintaining the high volt- 
age. Section 4.6 evaluates possible sputtering damage to the 
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large areas of the membranes and coatings, 
conclusions in Section 4.7. 


We present our 


4.2 DESIGN SUMMARY 

The ECCM concept has been described elsewhere by 


Mihora and Redmond. 


1151 


The basic strategy is to use elec- 


trostatic pressure 

t_ - e- l(E»n) E - 


E^n/21 


r 


f 


to form and (through electronic feedback) finely adjust the 
shape of a conductively-coated membrane. The nominal oper- 
ating conditions are a 50 kV potential across a 10 cm gapr 

2 

thus giving a pressure of 'vl nt/m , which is five orders of 
magnitude smaller than atmospheric pressure, but five orders 
larger than solar radiation pressure. The baseline system 
is a 100 m diameter antenna with positively biased, concen- 
tric ring electrodes. 

For future reference we set down here, in simplest 
form, some relations between membrane tension, shape, and 
applied pressure. The vector diagram below indicates the 
forces on a small circular piece of membrane. The electro- 
static force, TTr^tg.n, is opposed by the upward component of 
the membrane tension t^^^ (Figure 4.1) . Thus, if the radius 


f 


t 


Figure 4.1. 



Force balance vector diagram for a small section 
of membrane . 
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of curvature of the membrane is R . we have 

w 

irr^ . h ■ 2TTrt_ ain6 z 2rrt_ 

• lu in 

C 




■/V 2 

For the system under consideration, t.*n » 2Nt/m , R_ ~ 100 m^ 

'V C w 

so the membrane tension is t„ z 100 nt/m. 

Some drawings of possible ECCM designs are shown in 
Figures 4.2 and 4.3. All have the same basic features: 

e The high**field region is bounded by the membrane 
above, the electrode below, and some circumferential mounting 
structure . 

• The region behind the high-voltage electrods) con!- 
tains dielectric support structures. There may be substan- 
tial electric field in this region as well. The rear of 
the structure is also covered with membrane material. 

• An electronics package above the center of the 
antenna is held by booms mounted on the circumferential 
mounting structure. 

• Solar panels and other spacecraft hardware are 
mounted behind the antenna structure. 

The emphasis of this analysis' is on the high-field 
region and on the membrane separating it from the plasma 
environment* Also of primary importance is the amount of 
plasma leakage near the membrane circumference. 

The candidate membrane materials are 0.0005 inch 
(13 Jim) mylar with both sides aluminized to 1 fi/O surface 
resistivity (about 350 A of aluminum) r 0.0003 inch (8 ym) 
Kapton (presently available with one side aluminized) ; and 
various forms of metallized cloth. The material thicknesses 
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are such as to respond to the low values o£ pressure pro- 
vided by the electric field without being excessively dif- 
ficult to handle* 

Systems designers must take into consideration the 
thermal control aspects of having a large area of doubly 
aluminized membrane. While from an electrical point of 
view such material has great advantage in ensuring against 
loss of metallization, possible thermal problems may well 
need to take precedence. 

4.3 ENVIRONMENT SUMMARY 

We were asked to consider interactions of the elec- 
trostatic antenna with the space environment in both low- 
earth orbit (LEO, 'V200 km altitude) and geosynchronous 
orbit (GEO, 'V'? R altitude) . We did not consider the pos- 
sibility of polar orbit. While the ambient plasma in space 
is very tenuous, it is not negligible. Even in GEO, where 
the density is about 1 cm”^ , environment interactions have 

caused spacecraft malfunctions and even permanent damage. 

4 6 

In low-earth orbxt the plasma is colder and denser (10 -10 
—3 

cm ) and has even greater potential for causing problems 
to high-voltage systems. 

In Tables 4.1 and 4.2 we indicate the range of para- 
meters characterizing the low-earth and geosynchronous 
environments. While the SCATHA spacecraft is returning 
additional data on the GEO plasma, it is already clear 
that this environment is extremely variable and any design 
must take its full range into account. In interpreting 
Table 4.1, it should also be kept in mind that characterizing 
the plasma by a temperature tends to underestimate the high 
energy portion of the spectrum. For a thorough review, see 
Reference 16. 
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Table 4.1. Geosynchronous 

Plasma Environment 


Typical 


High 


5 -S 

Plasma Density 10 m 

Electron Temperature <1 keV 


Ion Temperature 
Electron Current 
Ion Current 
Debye Length 


1 keV 
lO""^ A/m^ 
10 "® A/xt? 


10® m"® 

, 1 keV 
5 keV 
10 “® h/'o? 
3xl0"® A/xt? 


5x10® m"® 

10 keV 
20 keV 
3xlo"® A/xx? 
lO"^ A/m^ 
10® m 


Figures exclude particles 

with e < 50 eV. 

r^— ■ .. - — 1 

Table 4.2. 

Low-Earth Orbit 


Environment 

Plasma Density 

10^® - lO®"® m"® 

Neutral Density 

- 10l« m-3 

Plasma Temperature 

~ 0.1 eV 

Debye Length 

~ 1 cm 

Spacecraft Velocity 

7500 in/sec 

Ram Ion Current 

- lO"® A/m® 

Ram Ion Energy 

- 5 eV 









The LEO (up to 1000 km) plasma has been studied for 
many years by electromagnetic, rocket and satellite probes. 

A thorough summary is available in Reference 17. The plasma 
density has a maximum of 10^-10 ■ cm in the range of 300 
to 400 km. It drops off sharply below 250 km. Above 500 km, 
it varies slowly with altitude, maintaining a value of 
10^-10^ cm”^. The plasma temperature takes a value of 
'^^0.1 eV from 200 to 500 km, though some measurements are 
as high as 0.3 eV. Below 200 km, the temperature decreases 
to '>^0.05 eV at 100 km. 

An additional feature in LEO is the presence of ram 
ions. Since the spacecraft velocity ('^'7.5 km/sec ) exceeds 
the ion velocity ('vl km/sec) , the forward side of the space- 
craft is constantly bombarded with '^5 eV ions and neutrals, 
which may, in some instances, cause sputtering problems. 

Photoemission due to sunlight is another environ- 
mental factor to be considered in both LEO and GEO. Most 

-5 2 

materials photoemit '\j2x 10 A/m of low energy (*^2 eV) 
electrons when exposed to solar ultraviolet radiation. 
Photocurrent is particularly important in GEO, where it is 
usually more than sufficient to maintain the spacecraft at 
positive potential. The excess photoelectrons then form a 
sheath of relatively dense plasma around the spacecraft, 
capable of supporting fairly considerable parasitic currents. 

In LEO, however, the photocurrent is considerably 
smaller than the electron thermal current and comparable to 
the ion (thermal or ram) current. Thus, the spacecraft will 
normally remain at a slight negative potential in sunlight, 
and no photoelectron sheath will form. However, when a 
spacecraft has applied differential potentials, photoemission 
should still be considered as a possible source of parasitic 
current. 
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4.4 EFFECTS OF SPACECRAFT CHARGING 

Spacecraft charging is an effect which occurs fairly 
frequently at geosynchronous altitude when a satellite en- 
counters a "geomagnetic substorm". During such events, a 
very hot (“VlO keV) plasjna can cause uniform (negative, rela- 
tive to plasma ground) or differential potentials comparable 
to the plasma temperature. The differential potentials can 
lead to arcing, causing malfunctions by, or damage to, space- 
craft. 

To prevent differential charging, the ECCM is designed 
to be conducting over its entire surface. This is accom- 
plished by covering the antenna structure with membrane mate- 
rial, which we will assume to be aluminized plastic. Thus, 
in analyzing spacecraft charging effects, we must concen- 
trate on what flaws are likely to develop in this scheme and 
what effects these may have on the antenna performance. 

The type of charging not eliminated by the "Faraday 
cage" structure is overall (uniform) charging. If the 
antenna charges uniformly to -10 kV, one would expect fields 
'^/10'^ V/m near the antenna center and '^^10^ V/m near its peri- 
phery. These fields will create small electrostatic pres- 
sures (<10 nt/m ) opposing the applied pressure ('\>1 nt/m ) . 
with the proposed feedback system, these stray pressures 
will be automatically compensated for by imperceptible 
(<0.1 percent) changes in the electrode voltages. 

Differential charging within the membrane due to 
buildup of trapped charge remains a possibility. If we 
consider a doubly aluminized plastic membrane of thickness 
10 ^ m, incident electrons with energies in the range of 
1 keV < e < 20 keV will be trapped. If the current of such 

~ ~ _5 2 

electrons is 10 A/m , internal potentials will build up 
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at the rate of a few volts/sec , or internal electric fields 

5 I 

at a few times 10 volts/cm-sec. Since the intrinsic 

-18 I 

conductivity can be 10 mho/m or less, electric fields of I 

8 i 

10 V/m may be required to balance the incident current. 

Such fields are in excess of the material's dielectric 
strength. High-field and radiation-induced conductivity 
will serve to reduce the likelihood of achieving breakdown 
potentials, but the possibility of occasional small arcs 
within the membrane must be acknowledged. These arcs will 
cause some material damage and some loss of metallization. 

Most likely f however, these effects will be sufficiently 
rare and sufficiently localized that they will not be signif- 
icant. 

Granting the possibility of some loss of metallization, 
we may address the question of how big a damage spot will 
degrade the antenna performance. A damage spot will have a 
substantial pressure reduction, and a displacement given by 
the square of the spot dimension divided by the radius of 
curvature: 

6 a 

c 

Thus, a 10 m displacement in an antenna with 100 m radius 
of curvature will result from a 10 cm damage spot. We 
believe development of such a large spot to be rather unlikely. 

Another easily avoided charging problem is differential 
charging of booms. Normal graphite-epoxy booms are suffi- ' 
ciently conductive to avoid charge buildup (required is 
surface resistivity less than 10 £?/□ ) while being resistive 

enough not to interfere with antenna function. Care must 
be taken, however, to avoid isolating a boom through use of 
insulating spacers or mountings. 



4,5 PARASITIC POWER LOSS 

An important neglected area that was not addressed 
in the preliminary design of the ECCM is the question of 
power loss. Any plasma entering the high-field region of 
the antenna constitutes a parasitic current and, when mul- 
tiplied by the applied voltage, a power loss* We have had 
little input as to what would be an acceptable power loss. 
Consequently, in this section we adopt the point of view 
that a kilowatt power loss is unacceptable, a one hundred 
watt power loss is a significant design consideration, and 
a power loss of ten watts or less is acceptable* 

There are three main avenues through which plasma 
(i.e., electrons) can enter the high-field region: 

e Through holes or gaps in the "Faraday cage.” We 
anticipate such gaps at the circumference of the membrane, 
where it is attached to the structure. 

e If metallized cloth is used, stray fields pene- 
trating the weave will enhance plasma collection through 
the membrane. 

• If a thin plastic membrane is used, high energy 
electrons will penetrate through the membrane into the 
electrode region. 

We consider the lattar two mechanisms in some detail below. 
First, however, we discuss the magnitude of current provided 
by the spacecraft environment. 

4.5,1 Environmental Sources of Current 

The first-order estimate of environmental current 
is the one-sided plasma thermal current: 
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Values for several environments are given in Table 4.3. Of 
these, the first three (low-density) cases represent quiet, 
moderate, and fairly severe geosynchronous conditions, while 
the two higher density cases bracket the LEO environment. 
Also shown is the typical photoemission current due to solar 
ultraviolet. 

In Table 4.4 these densities are multiplied 

by 50 kV to yield a power loss. The values in the electron 

(e~) column may be taken as reasonable upper bounds. From 

2 

these figures, it is apparent that in LEO even 1 m total 

4 2 

opening in a membrane of 10 m area is likely to produce 
an intolerable power loss. In GEO, however, a few square 
meters of opening may well be tolerable, provided enhanced 
collection of photoelectrons by stray fields is not too 
severe. (See Appendix C. ) 

Another possibility suggested by these tables is to 
bias the electrode negative. The ion currents which may be 
collected are far lower than the electron currents, and 
leakage of negative field would not lead to photoelectron 
Collection. In GEO, this would allow a great deal of lati- 
tude for openings and gaps in the membrane. It is unlikely 
that collected ions would cause sputtering problems. While 
negative bias may also be useful in LEO, the possibility of 
sputter damage for long-term missions must be considered. 

4.5.2 Field and Plasma Penetration Through Cloth Weave 

It has been suggested that metallized cloth be con- 
sidered as an alternative to aluminized plastic for use as 
membrane material. Metiailized cloth has been flown success- 
fully in the past; for example, the ATS-5 antenna was made 
of 70 percent transparent copper-clad Dacron. Suitability 


Table 4.3. Plasma Thermal Current 

Jq ■ eIk0/2 mj^ (A/m^) 


”o 


0 

e" 



1 cm"^ 

1 

eV 

2.7x10"® 

6.3x10"^® 

1.2x10"^® 

1 cm”^ 

1 

KeV 

8.5x10"^ 

2.0X10"® 

3.8x10"® 

1 cm"^ 

10 

KeV 

2.7xl0“® 

6.3x10"® 

1.2x10"® 

10^ cm"^ 

1 

eV 

2.7x10"^ 

6.3X10"® 

4xlo"® (ram) 

10® cm"^ 

1 

eV 

2,7x10“^ 

6.3X10"^ 

4x10“^ (ram) 

(SUNLIGHT) 


3x10"® 
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Table 4.4 . Power Loss (W/m ) 
0 50 KV 


"o 


0 


H^ 

4 

1 cm"^ 

1 

eV 

1.3x10"® 

3.1x10"® 

6x10 

1 cm"® 

1 

KeV 

4.3x10”® 

1.0x10“® 

1.9x10 

1 cm"® 

10 

KeV 

0.13 

3.1x10"® 

6x10 

10^ cm"® 

1 

eV 

13 

0.31 

0.2 

10® cm"® 

1 

eV 

1300 

31 

20 

(SUNLIGHT) 


~1 
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for rf applications does not, however, imply adeguate block- 
age of very high dc fields or impinging electrons. For this 
reason, we used a modified version of the fine resolution 
LEO current collection code to study field penetration and 
leakage of low-energy plasma through such material. 

Representation of two different grades of cloth mate- 
rial is illustrated in Figure 4.4. (The illustrated struc- 
tures are periodically contained so as to make a full mem- 
brane.) The two cloths have square cross-section fibers, 
with spacing-'to-diameter ratios of 8:1 and 4:1, giving 
nominal transparencies of 77 percent and 56 percent. We 
first calculated the potential (Figure 4.5) above and below 
the cloth in a configuration of high electric field below 
the grounded membrane, and a ground plane (possibly repre- 
senting a sheath edge) several fiber diameters above. We 
found the field above the coarser membrane 'v4 percent of 
that below, and correspondingly 1 percent for the finer cloth. 
From this we may expect the mean potential of a cloth membrane 
to be 'vl percent of the high voltage electrode, i.e. , 

^500 volts . 

We then proceeded to calculate the apparent transpar- 
ency of the membrane for low-energy electrons emitted from 
the ground plane. As seen in the particle trajectories. 

Figure 4.6, a substantial amount of focusing took place. In 
both cases the penetration exceeded the nominal transparency 
by about 20 percent. These results are summarized in Table 4.5 

These calculations lead to the conclusion that a sub- 
stantial fraction of electrons incident on a metallized cloth 
membrane will appear as parasitic current. Thus, there will 
be an intolerable power drain of several kilowatts due to 
photoelectrons alone. We are forced to conclude that metal- 
lized fabric is not suitable for use in the ECGM. 
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4.5.3 Plasina Penetration Through the Plastic Membranes 


The very thin (0.0003 to 0.0005 inch) membranes pro- 

k. 

posed for the ECCM will be penetrated by the hi<^h energy 
(c > 20 keV) electrons in the CEO plasma. While this sounds 
fairly high as a plasma temperature » half the current in a 
^ Maxwellian plas.aa has energy exceeding 1.7 times the tempera- 

ture. Purthermore? a thermal characterization of the GEO 
plasma tends to •'•nderestimate the high energy plasma component. 

Detailed theory and results for the power loss due to 

W this mechanism are given in Appendix P, We find that, for a 

-3 

plasma density of 1 cm , the power loss reaches 100 watts 
at a temperature of wio keV for the 0.0003 inch membrane. 

Thus, the system design for an ECCM in GEO must allow for 
• the fact that occasionally (<10 percent of the time) the 

antenna will require power in excess of 100 watts. A pos- 
sible solution is to shut down the antenrta when such a con- 
dition occurs. Such a solution, though painful, may be 
® preferable to maintaining the surplus power needed to assure 

continuous operability. 



63 



This type of power loss can, in fact, be reduced by 
spacecraft charging. Due to the conducting design, it is un- 
likely that the structure will ever charge to negative poten- 
tial in sunlight, since such potentials occur only as a conse- 
quence of differential charging. In eclipse, however, 
negative potentials will be attained in the same einvironments 
which produce high power loss. We estimate that an aluminum 
surface will charge to -10 kV, and a kapton surface to 
-9 kV in a 10 keV environment. This will reduce the power 
loss by nearly a factor of three, and should reduce the fre- 
quency of shutdown to a tolerable level. 

An interesting question is whether there are enough 

membrane-penetrating electrons in the LEO environment to 

cause substantial power loss. Even in the radiation belts 

(well above LEO) the current of high energy (>2 keV) elec- 

-7 2 

trons is less than 10 A/m , which would cause a power loss 
of only 'v50 W. Thus power loss due to high-energy electrons 
in LEO would be at most a few watts. 

An additional possible: power loss problem for these 
plastic membranes is back-side photoemission by transmitted 
ultraviolet. While we suspect that such photoemission is 
negligible f experimental verification would be worthwhile. 





4.6 


SPUTTERING 


At altitudes of 200-300 km the velocity of a satel- 
lite in circular orbit is about 8x10^ cm/sec. In the sat- 
ellite freune of reference, ambient atomic oxygen and molec- 
ular nitrogen (or O"^ and N^) have kinetic energies of about 
5 eV and 9 eV, respectively. Such energies must be only 
slightly above, if not below, the threshold for sputtering. 
Near threshold sputtering yields are not well known, but 

available measurements of sputtering show yields less than 

_2 

10 at energies less than 1.5 times threshold. 

5 3 

Ion concentrations are of order 10 /cm and neutral 
oxygen concentrations are about 10^-10^® cm~^. Assuming a 
yield of 10 , the time required to sputter a monolayer 

from a surface by species of concentration N (cm~^) is 


10 


15 


^ 

® 8x10® X 10"^N 


sec 


10® sec (ions) 

9-3 

100 sec (atomic oxygen at 10 cm ) . 

Thus, since the yield is not well known, sputtering of 
external surfaces by atomic oxygen in LEO poses a problem 
of uncertain proportions. Unless pertinent data on the 
sputtering of surfaces by 5 eV 0 atoms can be found, defin- 
itive experiments should be performed to assess the magnitude 
of the problem. 

Ambient ion fluxes of about 10^® cm”^ sec”^ may reach 
a negatively biased electrode through any openings in the 
Faraday cage. Assuming that the sputtering yield for 50 KV 
ions is about 1, the time for monolayer removal is about 

5 

10 sec., or 1 day. Thus, one year or so would be requxred 
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to scrub 1000 A of exposed areas of electrode surfaces, but 
the areas involved would be limited to about the size of the 
openings in the Faraday cage. 

4.7 CONCLUSIONS 

We conclude that, insofar as spacecraft-environment 
interactions are concerned, the electrostatically controlled 
membrane mirror (ECCM) is a viable concept for space applica- 
tions. However/ great care must be taken to enclose the high 
voltage electrodes in a Faraday cage structure to separate 
the high-voltage region from the ambient plasma. For this 
reason, metallized cloth is not acceptable as a membrane 
material . 

Conventional spacecraft charging at geosynchronous 
orbit should not be a problem provided ancillary structures 
(such as booms) are given non-negligible conductivity and 
adequate grounding. Arcing in the plastic membrane will 
probably occur occasionally for one-sided aluminization, and 
rarely for two-sided aluminization. However, loss of metal- 
lization sufficient to degrade antenna performance is highly 
unlikely. 

Power loss due to plasma electrons entering the high- 
field region is a potentially serious problem. In low-earth 
orbit any opening whatever in the Faraday cage is likely to 
produce an unacceptable power drain. At geosynchronous alti- 
tude, where current levels are lower, a gap of a,! cm at the 
membrane point-of-attachment may be tolerable. However, a 
geosynchronous antenna will occasionally encounter high 
energy plasma capable of penetrating a thin membrane to pro- 
duce power drains of one hundred watts or more. Provision 
must be made for shutting down the antenna during such 
periods. 
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The most credible sputtering problem is degradation 
of aluminized surfaces by ram neutrals in low-earth orbit. 
Even very low sputtering yields (<10 ) can seriously damage 

these surfaces on a time scale of days or months. If perti- 
nent data is not available, definitive experiments should be 
performed . 
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0 . B> P«rks^ M. J. Mandall and I. Katz 
Syitama, Seianea and Seftwara 
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Abitract 

Tha pr.rposa of tha praaant study is to 
datasmina tha capability of a fluid modal 
of alactron transport to axplain obsarvad 
propartias of ion thrustar ganaratad plas- 
mas. Calculations raportad hara show that 
whan tha affaetiva collision fraquancy in 
such a modal is Of tha ordar of tha alac- 
tron plasma fraquancy, tha rasulting alac- 
tric potential variations and alactron 
tamparaturas are in qualitative agreamant 
with values maasurad in tha plasma gene- 
rated by thi SEAT ZZ thrustar. Both theory 
and probe maaauramants made in flight and 
ground teats indicate substantial dapar- 
turas from tha barometric law and strong 
variations of plasma potential across the 
beam boundary . 

Nomanclatura 

B magnati,; field 

C collision operator in Boltzmann- 

Vlasov aquation 
E ^ ^ olaotrio field 
f drdv number of particles of type j in 
^ phaja^spaca volume element dr dv 

^ at r,v 

F total agergy flux 

^a ■ *^a^® 

I unit tensor 

} net currant density 

!c Boltzmann ' s constant 

m. mass of particle of species a 

a 

m electron mass ^ ^ , 

n electron density » /f(r,v,t)dv 

p ■ 1/3 n ft v'^dv ■ scalar electron 

, pressure 

P m ft v'v'dv ■ pi + 1, pressure 

tensor , 

Q ■ 1/2 mfv'^C dv 

q, charge on particle of species a 

A 

q magnitude df electron charge 

^ heat flux, Eq. (7) 

t position vector of a particle 

H * /mv C dv 

J electron temperature 

Xelosity of a particle 

Y" » V - V 

V mean ^or^dgift velocity ■ 
l/n /f (r,v,t)dv 

X mean free path for pair collisions 

^ between electrons 

0 » kT 

n plasma resistivity 

< thermal conductivity of plasma 

V effective collision frequency 
electron-ion collision frequency 

tr » nm W - <V >/3 I ■ stress tenson 


CofjfHflil ^ AmtHeflii Insiitvle of Acrontiilics mi 
Asirofimiks, (nr.. I9SI. All Hf dis mtrvti. 


electric potential 
electron plasme frequency 


1 . Introduction 

The purpose of the present study is to 
determine the capability of a fluid model 
of electron transport to explain observed 
properties of ion thruster generated plas- 
mas. Calculations reported here show theit 
when the effective collision frequency in 
such a model is of thii order of the elec-^ 
tron plasma frequency, the resulting elec- 
tric potential variations and electron 
temperatures are in qualitative agreement 
with values measured in the plasma gener- 
ated by the SEHT 12 thruater. Probe mea- 
surei^en^s made in $ERT XI flight and 
ground test experiments (2) indicate sub- 
stantiai departures from the barometric 
lawn~7) 2 ind show strong variations of 
plasma potential across the beam boundary. 

We propose to explain the plasma proper- 
ties observed in the aforementioned experi- 
ments in terms of anomalous resistance of 
the plasma to the flow of electron current. 
The calculations are based on fluid equa- 
tions expressing conservation of charge, 
momentum, and energy. We adopt the clas- 
sical (ignoring thermoelectric effects) 
form of. the equations of electron trans«* 
por^,'^^ but permit reduced values of the 
transport coefficients. Predicted space 
dependent potentials, electron tempera- 
tures and current densities agree qualita- 
tively with experimental results. 

While the plasma is not collision domi- 
nated, randomization of electron velocities 
may still occur through enhanced levels of 
fluctuating fields, such as those initiated 
by streaming instabilities. Such fields 
are probably effective in coupling neutral- 
izer electrons into the bulk plasma and in 
equalizing the mean drift of electrons with 
ions in the thruister beam. Such mechanisms 
are often approximated by introducing an 
effective collision frequency, v. 

Within a meter or so of the ion thruster, 
the electron densities are in the range 

10® < n < 10^^ cm**® 

and their velocity distribution is charac- 
terized by temperature 0 between about one 
and ten -slectron volts. The Debye length 

^ ZOoyT cm 

is typically small compared to distances L 
over which there is a substantial variation 
of macroscopic plasma properties such as 
density, potential, and temperatures. On 
the other hand, the mean free path for pair 
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collisions Xq 

X_ 3 10^^ 0^^^/n cm, E < 0 

C 

for alactrons of energy E (eV) is typically 
long compared to hf no that as previously 
asaerted the behavior of the plasma is con- 
trolled by collective rather than colli- 
sional effecti. Since << L, the plasma 
ii quaeineutrali departures from neutrality 
amounting roughly to 



the space around the vehicle is strongly 
shielded from surface potentials. This is 
in contrast to the situation that prevails 
in charging of spacecraft in geosynchronous 
orbit where effects of space charge are 
entirely negligible and potentials are 
determined as solutions of Laplace ^s equa- 
tion* 

Although collisionlessr thruster-gene-^ 
rated plasmas exhibit macroscopic behavior 
similar in many respecta to that of a col- 
lisional plasma. Such behavior is perhaps 
not totally unexpected in view of the fact 
that in both non-equilibrium and equili- 
brium plasmas electrons are scattered by 
fluctuating electric fields, A primary 
difference between the equilibrium and non- 
equilibrium cases is in the magnitude of 
the fluctuating fields. 

Several investigators have measured 
properties of thruster generated plas- 
mas. V3-7) the experiments of Ogawa^ 

et al . , on cesium ion beams neutralized 
by electrons from a hot wire, measurements 
were made of the density, potential, and 
electron temperature in the beam plasma. 

The potential difference between the 
neutralizer wire and the plasma could be 
varied by changing the position of the 
wire, the large potential differences 
(electron injection voltages) occurring 
when the wire was completely withdrawn 
from the beam plasma. An important result 
of the Ogawa experiments was that over a 
wide range of conditions electron density 
n and plasma potential were well corre- 
lated by the barometric law 

n(r) * const exp(q<|> (r)/kT) (1) 

The approximate validity of the barometric 
law was further verified by Kaufman.''' 

Since the barometric law is a thermal 
equilibrium concept, it can be completely 
valid only if the plasma is isothermal. 

The plasma is only approximately isotherm- 
al, noticeable deviation occurring as one 
proceeds from the beam axis beyond the 
beam: edge into the plasma formed by am- 
bient and charge exchange ions. Kaufman 
observes an electron temperature in the 
charge exchange plasma only |bout half 
that in the beam . ' ' ' Ogawa f and 
Sellen'3) obtained measurable temperature 


variations in the beam plasma over several 
tens of centimeters in the downstream direc- 
tion from the accelerator grid. The largest 
deviations from the barometric law were 
observed for large injection potentials 
(^10 volts) . Probe traces in such cases 
also indicated departures of the electron 
spectrum from a Maxwellian shape. 

Probe measurements of the plasma poten- 
tial in the thruster beam were made in 
SERT IX flight and ground test experiments. 
The measurements show strong variation of 
plasma potential across the beam boundary 
about 20 cm downstream from the thruster 
grids. Such results are difficult to ex- 
plain on the basis of a barometric law 
relationship unless the electron tempera- 
ture or density variation from beam center 
to beam edge is much higher than might be 
expected from other measurements made in 
similar configurations. We anticipate, 
however, that such is not the case and, in- 
stead, that the observed behavior should 
be explained in terms of the anomolous 
resistivity of the thruster generated plas- 
ma to the flow of electron current. Thus, 
the primary objective of the following 
sections of the report is to determine the 
capability of simple transport models to 
explain, at least qualitatively, the ex^ 
perimental results. 

The next section summatizes the kinetic 
equation for the electron distribution and 
the first few moment equations expressing 
conservation of charge, momentum, and 
energy. In Section 3 we state the approxi- 
mations leading to the transport equations 
which we eventually solve, The method of 
solution and the results of calculations 
will be the subject of Sections 4 and 5 
respectively. The final section, Section 6, 
summarizes the conclusions of this study. 

2. Exact Equations for Electron Gas 

In principle, a kinetic approach based 
on the Vlasov-Boltzmann equation fully 
describes the spacecraft generated plasma. 
The complexity of such an approach, how- 
ever, makes it impractical as a basis for 
conducting multidimensional calculations 
of plasma behavior. Besides, except near 
sources and collecting surfaces, where the 
distribution function may change markedly, 
one should be able to adequately describe 
the plasma in terms of certain average 
properties of the distribution, such as 
temperature, density, and particle and heat 
fluxes. Below, the exact equations de- 
scribing the plasma are given in order that 
the reader may be aware Of the effects 
neglected in arriving at the approximate 
equations that are subsequently solved. 

Quite generally the state of the plasma 
can be sgegified by the distribution func- 
tion fa(r,v,t) that characj^ejize gach parti- 
cle component a, where fa{r,v,t)dr dv 
represents the number of particles of 
species a^in^the six dimensional yo^ume 
element dr dv about the position r,v in 
phase space. The kinetic equations which 
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d«ierlb« th« diitributlon «r« 

it. - K 

m * V*7f ¥ • 7 f M C 

y— -r V t y^ 


(?) 


For particle* of m§ea m« and ehec?* in 

i n electric field E end e megnebic field 
the ** smoothed" force on e perticle is 
F«. The effecta of collisions between 
perticle* is taken account of by the col- 
lision term denoted her* by C^. Her* w* 
attempt to describe the plasmi in terms of 
its density n, mean velocity V, and cer- 
tain higher velocity moments. For conven- 
ience, w* have omitted the particle 
species subscript a. The first three 
moments of the kinetic equation yield 
conservation equations for particles, 
momentum and energy, as summarised be- 
low»(l) 

Conservation of Particles 

7*nV • 0 (3) 

0 w 

ConiTvation of Momentum 

(mn^) + 7*mn^ + 7*P - qn/f + 

' (4) 

Coniervation of Energy 

■ qnE*V 

+ R*V + Q (5) 


electrons as a whole and that tlnsir distri- 
bution varies Slowly in space* 

3 . ADDroximations for Electron Oaa 

Consider that the plasma is in a steady 
stats end that quasi-neutrality pertains 
throughout the bulk plasm* (that is, away 
from electrodes and collecting surfaces) . 
The electrons and ions each satisfy the 
particle continuity equation 

7*n^V^ ■ 0 (i ■ +*") (9) 

with 04 - • n. ■ n. The momentum equation 
simplifies considerably if the electron 
drift velocity V is small compared to the 
random velocity <v'"*>^/2 and if the veloc- 
ity distribution is nearly isotropic. 

Then, in the absence of magnetic fields, 

7p + eni ■ 15 (lO) 

where R represents the collisional drag be- 
tween ions and electrons. In a glassical 
plasma dominated by colliaions, R is com- 
posed of a p§rt proposticnal to the rsla- 
tive motion u ■ V« - Vi between electrons 
and ions, leading to plasma resistivity, 
and to a thermal part proportional to the 
gradient of electron temperature, which is 
frequently neglected. In this approxima- 
tion, equation (4) becomes 

7p + enE ■ nne^ (11) 

where ? is the net current density and the 
plasma resistivity n is related to the 
electron-ion collision frequency v«i by 


where 
F ■ 

is the 

q - 

is th« 



1 

jjv + ifV -f q 

-1 “o 1 

L 

total energy flux* 

•* (6) 

If the plasma is non-resistive and iso- 
thermal* equation (11) yields the baro- 
metric law* equation (ij. In this sense* 
equation (11) * or more generally the com- 


(7) 

plete electron momentum equation* may be 
regarded as the generalization of the baro- 

heat flux* 


metric law. 


Q 



(B) 


and < > denotes an average over the dis- 
tribution f. 


So far* the equations are quite general 
and involve no assumption that the gas is 
collision dominated or retains a Maxwellian 
spectrum of velocities. Separate conserva- 
tion equations may be written not only for 
different particle species ^ but also for 
different groups or particles of the sam^ 
charge and mass. Primary electrons* for 
example, with significant streaming ener- 
gies could be treated as distinct from the 
main electron population which is taken to 
have a Maxwellian distribution of veloci- 
ties. For the present* however* and until 
experimental or theoretical considerations 
dictate otherwise* we shall consider 


If the plasma is not collision dominated* 
randomization of electron velocities may 
still occur through the enhanced levels of 
fluctuating fields in the plasma* such as 
occur for electron two-stream instabili- 
ties* or electron-ion instabilities of the 
ion-acoustic or Bunemann type.'®/B^ These 
mechanisms are probably effective in coup- 
ling neutralizer electrons into the bulk 
plasma and in equalizing electron and ion 
mean drift velocities. They are often 
approximated by introducing an effective 
collision frequency* v* in place of 

The determination of electron tempera- 
tures in the plasma requires considera- 
tion of the energy balance equation* equa- 
tion (5). Making the same approximations 
in the equation expressing conservation of 
energy that were made in the momentum 
equation* yields 


3 
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(13) 


V*? • qnE*V + 5**^, + 

with 

F - I + q (14) 

Here 5/2 it the tQthalpy flux of th« 
drifting •Itetronj ^ q tht macrotcopic 
h«at flux, and R»V ia ralatad to tha af- 
factiva joula haating aaaociatad with tha 
relativa motion of alactrona and iona. 

Tha quantity R appaara alao in tha alac- 
tron momamtuffl aquation; for a plaama con* 
trolled by collactiva affacte it ahould 
be approximatad in tha energy aquation 
in the aame manner a a in t||{e momentum 
equation* Tha heat flux, q, containa new 
faaturea. Claaaically, q containa two 
terma; one proportional to tha relative 
drift velocity between alactrona and iona, 
and tha other proportional to tha gradient 
of electron temparatura. (8) 

For the initial calculationa, we ignore 
the drift contributiona to the energy flux, 
the electron^ion haating Qei, and aaauma 
that tha heat flux ia proportional to tha 
temperature gradient. The energy balance 
equation thua aaaumea tha aimpla form 

V»k 70 + mnv * 0 (15) 

4. Ion Engine Neutralizer Code 

The baaic physics of an ion engine 
neutralizer model was presented in the 
preceding sections. This physical model 
has been incorporated into a two-dimension- 
al (R-Z) computer code, which is described 
below. A sample calculation of neutrali- 
zation in a thruster similar to SERT IZ 
is discussed below. Results are given 
for space dependent electric potentials, 
electron temperatures, and current densi- 
ties. 

Code Description 

The ion thruster model has been incor- 
porated into a two-dimensional (R-Z) com- 
puter code following the block diagram 
shown in Figure 1. For this initial 
version the ion currents and densities 
were assumed known, in a later version, 
it would be possible to allow a multi- 
component ion composition to be determined 
self-consistently with the temperature and 
potential. The code operates entirely in 
MKS units. 

The code has been tun interactively, 
with all relevant information on disk file. 
As long as previous information exists on 
disk, the program may be entered from the 
two noted entry points as well as the be- 
ginning, For developmental purposes, it 
was found convenient to "hard-wire" many 
features of a particular problem into the 
code, while others are prompted for input. 

A flexible on-line graphics program, which 
plots information on the disk file, has 
also been developed. 



Fig. 1 Block diagram for ion engine 
neutralization code. 

At present, the code assumes ion veloc- 
ities everywhere and ion currents at the 
input boundary to be known. The code then 
calculates plasma densities such that 
7*(ny) ■ 0 is numerically satisfied. Typi- 
cally, ion velocities are taken to be either 
purely axial or to be radial from a point 
source on the axis ekterior to the mesh. 

As the code requires non-zero plasma den- 
sity everywhere, a background density of 
"slow" ions may be added. It should be 
possible to handle multiple ion species 
with interconversion fairly easily. 

The neutralizer is assumed to be a ring 
at specified distance from the axis, emit- 
ting a current of electrons equal to the 
ion beam current. The net current in the 
plasma is given by 

j ■ nq(v^ - Vg) (16) 

. a (-7 $ +2 ; (17) 

^here now p ■ nJcT is the electron pressure. 
(For o-^oD and qkt-^e (constant) we find 
* 9in n.) The code determines electro- 
static potentials by solving * 0. (See 
Appendix As) It is necessary to iterate 
between this equation and the temperature 
equation (equation 15 ) r since the pressure 
is a function of temperature. 

The plasma temperature satisfies the 
equation 
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a 


(18) 


wh«r« K 1$ th« 'ch«nn«l eenduetlvity and 
th« ri 9 ht->hand>tid« tapraianti th« ohmic 
9 «n*r«tion of h«at. For fehla pralimlnary 
varaioiif convactiva haab transport has 
baan naglactad. On tha various boundary 
ragions, althar isotharmal or Insulating 
boundary conditions may ba spacifiad. 

Sinea> in practiea# wa taka k to hava a 
poifar law dapandanca on df *c ■ nK'’9(*'**> 
tha aquation actually solvad is 

-7»k'7(6") - . (19) 


For convanianea« tha transport coaffi- 
ciants 9 and k" ara caiculatad by a singla 
isolatad subroutina. Tha conductivity a 
may dapand on both dansity and tamparatura> 
and k' on dansity only. Tha peasant var- 
Sion assumas a ralaxation rata proportional 
to tha plasma fraquaneyt 


** sr 



( 20 ) 


where the peremeter a ie taken to be 
QaSla By the claseical (^eideman*Franz 
lawp 


the bean from the aide. The ion current ie 
SO percent neutraliied at MS cm downstream 
from the thrusters 



Fig. 2 Plaama density for SERT XZ ion 
thruster model • 



( 21 ) 


TOmaTUSCS (eV) 


Zf we meaaure temperature in eV| k ■ qp 
•o that • 3/4 a. 


5« Computational Resulti 


A calculation was performed for neu- 
tralization of a 0.23 ampere p purely axial 
beam of 3 keV mercury ions in a constant 
density background plasma. The beam had 
a radius of ^7 emp and a ring neutralizer 
was placed at a 17 cm radius. The given 
plasma density (Figure 2) had a peak of 
16.1 X 10^^ and an ambient density 

of 1.0 X loi^ m* These conditions ap- 
proximate those occurring in the plasma 
produced by the SERT XZ thruster/ The 
temperature profile (Figure 3) was calcu- 
lated with insulating boundary conditions 
at the thruster. The maximum temperature 
occurs at the beam entrance p where the 
heat generation is greatest. The peak 
temperature was 5 eVp compared with a 1 eV 
assumed background. The electrostatic 
potentials are shown in Figures 4 and 5. 
Strong potential variations across the 
beam edge and a potential dip near the 
neutralizer are calculated features in 


qualitative agreement with experimental 
results (Figure 5) . While the strong 
edge fields conform approximately to the 
barometric law at the local temperature p. 
they deviate substantially from the results 
that would be obtained with u << Up or by 
using the barometric law as the point of 
departure (that is, a zero resistivityp 
isothermal calculation) . Current-vector 
plots (Figures 6 and 7) indicate that the 
beam is neutralized by electrons entering 


Fig. 3 



Temperature profile for SERT II 
ion thruster model . 
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6. Summary and Conclusions 

The consequences of the assumption that 
the electron gas near an ion thruster be- 
haves as a resistive fluid have been ex- 
amined. Theoretical results obtained here 
are in qualitative agreement with experi- 
mental observations. Such agreement indi- 
cates that the properties of thruster 
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gcntrAtad plawMS cm b« d«fcrib«d by 
fluid uguationa having a elaaaieal form 
bu« with an affactiva eolliaion fraguaney 
naar u-, muoh in axcasa of tha olaaaioal 
valua xor pair (alaetron-ion) eelliaiona. 

Furthar work ahould ba parfomad to 
taat gualitativa and quantitativa pradie* 
tiva eapabilitiaa of fluid modala of 
thruatar plaaaaa« and to battar undaratand 
tha ralationahip batwaan ariata fluid 
modala of tha typa invokad hara and tha 
undarlying plaama phyaical maehaniama 
ambodiad in tha colliaionlaaa Vlaaov 
aquation. 

kpoandix A. Mathamatieal Conaidarationa 
on tha variational Formulation of 
boiaaon’a touation 


whara 

Wjij ***- 

d*r MjCr) #(r) , (A.<\ 

M 4 (f) ara nodal intarpolation funetiona 
and f i ara noda valuaa for tha unknown 
fiald* 


finea W ia ayimnatrie, thaaa aquationa 
ara ^ 

EHjj - 0 (A. 7) 

Uaing (A. 5) wa find 


Zn our thaory of ion angina nautraliaa* 
tion» it ia naeaaaary to aolva aquationa 
of tha form 



7-to(r)7 4) - s{r) (A.1) 

aubjaet to fixad-valua boundary eonditiona 
at soma nodaa and normal*gradiant boundary 
eonditiona at othara. ha naad to ahow 
that aquation (A.l) ia axaetly tha aqua* 
tion aquivalant to tha variational formula- 
tion » and that tha normal-gradiant bound- 
ary eonditiona ara aquivalant to a aurfaca 
eharga. 


Thaoram 1 

Minimization of 


/d^r 

0 ’ 


, .2 

-3^ l!^l 


+ <Mr)S(r) 


(A. 2) 


ia aquivalant to aquation (A.l) « 


Whara 


?(r) ■ 

tha intagrand of (A. 8) can ba rawrit- 

tan as 

7‘{c(r)H.{r>7?(r)} 

- N, (r)7.(0(r)7?(r)l (A. 9) 

Writing Q(r) ■ o(r)^7(r)f and using tha 
divarganea thaoram > aquation (A. 7) ba- 
oomas 

fd^r Ni (r) D(r)*n -/d^r Mj (r) 7*D(r) 

axg i ^ ^ •• * •• 


Proof 

Minimization o£ fd^tL(f,V^,r) raquiras 
! ‘ Ti^r * — 7 $"" ' " ' ■ ° 

Taking L as tha intagrand of (A. 2 ) » aqua- 
tion (A. 3) yialda axaetly aquation (A.l). 

Thaoram 2 


+ - 0 (A. 10) 

whara tha first intagral runs ovar tha 
surfaca of voluma 0. By analogy with 
alaotrostaties. tha sacond and third 
tarma ara tha "voluma eharga" aasociatad 
with noda i, whila tha firsts involving 
tha surfaea-normal-gradiant of tha potan- 
tialf is non-zaro only on surfaoa nodaa 
and may ba companaatad by a eorrasponding 
surfaea-eharga density • 


In the finita-alamant formulation for 
the minimization Of (A. 2), normal -gradient 
boundary eonditiona ara aquivalant to sur- 
faea eharga. 


Proof 

Zn the finite-alamant formulation, 
minimization of (A. 2} leads to the aqua- 
tiona 


3 

3^1 



«jk ♦ 


E 

t 




0 

(A. 4) 


Corollary 

Consider tha aquation 

7*J ■ 0 (A. 11a) 

whara 

J ■ o(7<j>-g(r)7f) , (A. 11b) 

and f(r). g(r) are known functions. If 
the inhomoganaoua term is evaluated using 
aquation (A. 8), aquation (A. 11). with no 
further qualifications, will give zero- 
normal -cur rent boundary eonditiona when 
solved by the finite-alamant method. 
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Furthtrmor*, non-jstto^normil ourrtnt 
boundary conditional n 3^{x) , can ba in* 
vokad by adding an inhomoganaoua tarm to 
tha aurfaoa noda aquations » 

Proof 

Substituting (Aili) into (A<10) i tha 
finita aiamant aquations ara 

/d*z Ni(r) J (r)*n - /d%* j(f) Ni(r) - 0 

(Alia) 

Ab tha lacond tarm is to vanish (within 
tha finltaoalamant approximation ) , it fol- 
lowB that tha first tarm is aqual to tha 
right-hand^aida. ThuSf raplacing tha 
right-hand-iida with 

/ d^r M,(r> J.(r) (A, 13) 

will product tha apaoifiad boundary condi* 
tiona, 
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APPENDIX B 

MATHEMATICAL CONSIDERATIONS ON THE VARIATIONAL 
FORMULATION OP POISSON'S EQUATION 

In our theory of ion engine neutralization, it is neces- 
sary to solve equations of the form 




(B.l) 


subject to fixed- value boundary conditions at some nodes and 
normal-gradient boundary conditions at others, we need to 
show that Eqn. (B.l) is exactly the equation equivalent to the 
variational formulation, and that the normal-gradient boundary 
conditions are equivalent to a surface charge. 

Theorem 1 

Minimization of 

/d^r {2^1v<|)|2 + ,})(r)S(r)) (B.2) 


is equivalent to Eqn. (B.l). 

Proof 

Minimization of J*d^rL ((|) ,2'T?/r) requires 

n . _ 3 L (<{> ^ « (B . 3 ) 

- 3(2<j)) 34) ~ ^ 

Taking Las the integrand of (B,2), Eqn. (b. 3 ) yiei<is exactly 
Eqn. (B.l) . 
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Theorem 2 



In the finite-element formulation for the minimization 
of (B,2), normal -grad lent boundary conditions are equivalent 
to surface charge. 

Proof 


t 




€ 


i 


In the finite-element -formulation, minimization of (B.2) 
leads to the equations 


^<t>. 




lit "jk “ 


(B.4) 


where 


”jk ’ /n (e)l ■ 1! 11,^ (E)] (B.5) 


■ /n s (r) , 


(B.6) 


are nodal interpolation functions and 4>j_ are node values 
for the unknown field. 


Since W is symmetric, these equations are 


? "i: +3 * = » 


(B.7) 


Using (B.5:) we find 

Swij d^r [7 (r) i ' [0 (r) V? (r) ] 


(B.8) 


where 


(j, (r) = (r) tf)^ 
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The integrand of (B.8) can be rewritten as 

y • [a(r)N^ (r)y^(r) ] - (r) 7 • [a(r)y$'(c)] (b.9) 

Writing p(£) » CT(r)V^(r), and using the divergence theorem, 
Eqn, (B,7) becomes 

/g d\ N.(r) D(r) -ft - d^t N.(r) y-R(r) + - 0 

(B.10) 

where the first integral runs over the surface of volume 
By analogy with electrostatics, the second and third terms 
are the "volume charge" associated with node i, while the 
first, involving the surface-normal-gradient of the potential, 
is non-zero only on surface nodes and may be compensated by 
a corresponding surface-charge density. 


Corollary 

Consider the equation 

V • j = 0 (B.lla) 

where 

j = a(V(j)-g(r) 7f) , (B..llb) 


and f(r), g(c) are known functions. If the inhomogeneous term 
is evaluated using Eqn. (B.8), Eqn. (b. 11) , with no further 
qualifications, will give zero-normal-current boundary condi- 
tions when solved by the finite-elestent method. Furthermore, 
non-zero normal current boundary conditions, ft Jg(r), can be 
invoked by adding an inhomogeneous term to the surface node 
equations. 
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Proof 


Substituting (B. 11) into (b. 10) , the finite element 
equations are 

JV d^r N^(r) J(r) -ft - /d^r V-J(^) ) = 0 (B.12) 

As the second term is to vanish (within the finite-element 
approximation), it follows that the first term is equal to 
the right-hand-side. Thus, replacing the right-hand-side 
with 

/g d\ N. (r) J 3 (r) (B.13) 

will produce the specified boundary conditions. 


APPENDIX C 


ELECTRON PENETRATION OP THIN FILMS 

A be 2 un of electrons incident upon a solid surface will 
penetrate the material and \indergo collisions with both its 
nuclei and electrons. The nuclear collisions lead to large 
angle scattering and attenuation of the beam in the incident 
direction. The electron scattering causes less deflection ^ 
but reduces the translational energy of the individual elec- 
trons. Eventually, when the energy of the incident, or primary 
electrons has been reduced to (effectively) zero, they reach 
their maximum depth of penetration and stop. This maximum 
depth is called the "range.” The range, R, is often consid- 
ered to depend only upon the initial incident energy, E^, and 
may be represented: 

n, n- 

R a b, E^ + b, (1) 

1 O 2 0 

Some of the energy absorbed by the material causes 
secondary electrons to be excited into conduction bands and 
to migrate to the surface and escape. For a material thick 
compared to the range, the escape-surface is always the same 
as the incident surface. An empirical theory of secondary 
emission is well established for this case. 

The secondary electron yield, 6, (or ratio of secondary 
emitted to incident current) is given by the integral of the 
rate of production of secondaries at a. depth x , multiplied by 
their probability of escape to the surface, over all x. The 
empirical theory employs the following model: 



1. The rate o£ production o£ secondaries at a depth x 

is proportional to the rate o£ primary electron energy loss, 
or the "stopping-power" S(x), at that depth. 

S(x) * II . (2) 

2. The probability of escape of the secondaries once 
formed is proportional to e , where a is a constant, and 
1/a is a characteristic "depth-of-escape" for each material. 

The yield <S , for a normally incident beam is then 
given by; 

R 

(S(0) * C y* S(x) e“®^ dx (3) 

0 

For a beam incident at an angle 9, the distance to the sur- 
face is no longer x, but xcos9, i.e. 

R 

<S(0) * C f S(x) dx 

0 

Thin Films 

For a film with a thickness on the same order as, or 
less than, the range, both secondary and primary electrons 
will be emitted from the back surface as well as the incident 
surface. 

A central assumption in the empirical theory above 
is that secondary electrons are formed with isotropic initial 
velocities. Thus the probability of a secondary electron 
reachina the back surface has the seme form as that for the 
front siir race; i.e. , proportional to e (Figure 1) . 





Figure 1. 

The yield at the back surface is given by: 

R 

6^(6) * C y* S(x) dx (4) 

0 

The fraction of primary electrons that are emitted from the 
back surface n(9)f is given by the transmission coefficient. 
This depends on the nuclear scattering and the beam attenua- 
tion, discussed above. The fraction of the incident beam 

* "'NOX 

scattered in the forward direction is given by e , where 

0 is the nuclear total scattering cross-section, a depends 

on both the initial energy (and hence the range R — recall 

Eq. (1)— and the atomic number Z) . Using the parameterized 

theory of Kanaga and Okayama (J. Phys. D. 5 43, 1972). 

n(0) . (5) 
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For a film of thickness T, the path- length for a beam 
incident at an angle 6, through the film is T/cosd. 


n(e) . IT/(R COS0-T)] 


( 6 ) 


The total current emitted through the back surface is the 
sum of the transmitted primary and secondary contributions. 


W * ^in <’> 

or, writing the total yield, t(0) as: 

t(9) « 6^(9) + n(9) (8) 


i 


out 



T(9) . 


( 9 ) 


In a more realistic space environment the electrons are 
derived from an isotropic plasma, with a Maxwellian distri- 
bution of initial energies E^. In this case the total 
emitted current will be Eq. (9) averaged over all angles, 
and integrated over a Maxwellian distribution of initial 
energies. 

1 

t^®°(Eq) = 2 J" t(9) COS0 d(cos0) (10 

0 


. MAX 
out 



-B^T 


where T is the temperature of the Maxwellian and 


FO 



1/2 


where N is the density and m is the electron mass. 



(11 


(12 
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Power Loss for an Aluminized Kapton Film 


Calculation cf the emitted current from the back sur~ 
face of a film made of two materials A and B follows the 
above scheme very closely: 


Secondary Yield 


i. When the depth in the film D is less than the thick 
ness of the layer, A, presented first to the incident beam 
(outer) . 


S (x) 


dR^ (x) 


where x ■ D/cos6. 


6 (D) * S(x) e 


-la^(VD) + (.gTgl 


ii. When 0 exceeds secondary production occurs 
in the sedond layer and: 


S(x) » 


(15) 


where we calculate the remaining energy upon entering B(B^ ) 

OB 

via the range formula in A. 

A A A A 

■^1 ®0B *^2 ^2 ®0B “ \ ” (T^/cos0) (16 


6(D) » C« S(x) e 




To be transmitted the range R must exceed the total 
distance traveled inside the film, x » T/cos6, where T ■ + 

the total thickness of the film. The range, R, is calculated 
in two parts. 


i. Layer A 

If R^, calculated from (1) exceeds x^ * T^/cos0, the 
electron will proceed into layer B with energy Eqq given 
by (16) above. The fraction of the incident beam entering 
layer B is 

-0.187 <=os0 - T.)] 

- e ^ ^ ^ (18) 

ii . Layer B 

The range in layer B, R^, is given by the range 
formula using as the initial energy. The fraction 

transmitted rig is calculated via (18) above. 

The overall transmission coefficient is then 

n =» ri^ X Hg (19) 


Integrating over an isotropic Maxwellian gives the overall 
emitted current per unit area. For a 100 m diameter disk, 
above a +50 keV potential, the total power loss will be; 


I ISO 




50 X 10 W 


• w 




Results 


Using the above formulation the yields for a series of 
monoenergetic beams and the power losses for a series of 
Maxwellian environments were calculated for the following 
100 m diameter disks of films i 

1. 0,3 mil K75, 000 A) Kapton coated with ^<350 A of 
Aluminum, with the Kapton facing the incident beam. 

2. 0.5 mil ('\^125,000 A) Kapton coated with '^350 A of 
Aluminum, again with the Kapton on the incident sur- 
face. 

The yields and power losses are tabulated below. 
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POWER LOSSES 


0.3 mil (75,000 A) Kapton Backed with 350 & Aluminum 

Density ■ 1.0 cm“ 

Incident Current Emitted Current 
(u Am"^) (u Am”^) 


Temperature 

(kev) 


Power Loss 
(W) 



TRANSMZSSZOK YZELDS 


5 mil (125,000 A) 

Kapton Backed 

with 350 A 

Aluminum 

Bn«rgy (keV) 

Secondary 

Primary 

Total 

28 

0 

0 

0 

30 

0.023 

0 

0.023 

32 

0.053 

0.001 

0.054 

34 

0.063 

0.018 

0.080 

36 

0.064 

0.053 

0.12 

38 

0.065 

0.10 

0.17 

40 

0.064 

0.15 

0.22 

50 

0 . 053 

0.40 

0.46 

75 

0.035 

0.72 

0.76 

100 

0.026 

0.84 

0.87 
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POWER LOSSES 


0.5 mil (125/000 A) Kapton Backad with 350 A Altiminum 

Dansity ■ 1.0 cm"^ 

iMiENiratur. C»r«nt Emittid Current 

(kaV) (u Am"'^) (u Am (w) 


6 

2.1 

0 

0 

8 

2.4 

0.043 

17 

10 

2.7 

0.13 

52 

12 

2.9 

0.25 

100 

14 

3.2 

0.41 

160 

16 

3.4 

0.58 

230 

18 

3.6 

0.78 

300 

20 

3.8 

0.99 

•390 

30 

4.6 

2.0 

800 

50 

6 . 0 

3.9 

1500 
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POWER LOSSES FOR 1.0 cm DENSITY MAXWELLIANS 






APPENDIX 

Computer Code LEAK 


To calculate electron penetration and power loss through 
a general two-material filmr a small code (LEAK) has been 
written and has the following features: 


1. Subroutines 

a. LEAK - main program. 

b. TRASEC * calculates the transmitted secondary yield. 
TRASEC has the same argument list as ELSEC , etc . , 

in MATCHG and NASCAP. 

c. TRANIT - calculates the transmission coefficient, 
and has the seune format as TRASEC. 

d» DELC - calculates the integrand for the numerical 
integration in TRASEC. 

e. TCOEF - calculates the integrand for the angular 
integration in TRANIT, 

f. STOP - solves the stopping power equation (16) for 
Eqq and the stopping power. 

g. RANGE - calculates the range for each material given 
an initial energy. 

h. MAXWSR - integrates the yield over a Maxwellian. 
Identical to NASCAP and MATCHG routines. 


2. Input 

The present version is fixed to calculate yields only 
for isotropic incident electrons and power losses only 
for 100 m diameter disks above a +50 keV potential. A 
more flexible version is trivial to achieve if required. 

a. DATA INPUT - The four range parameters of each 
material and their secondary emission parameters 
(C, a) are supplied in a DATA statement. 
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b. INTERACTIVE INPUT 

- thickness of incident layer. 

- thickness of emitting layer, 

- temperature step size (T^-keV) . (20 values of output 

variables are determined at temperatures T, 2T... 

• 20T.) 

- the density of the Maxwellian (N^ m**3) . 

- the film potential in volts. 


3. Output 

Two tables are produced: 

a. A table of secondary/ primary^ and total yields for 
monoenergetic beams of electrons. The energies are 
given by T, 2T...20T. 

bu A table of incident and emitted current and power 
loss for Maxwellian environments with density N and 
temperatures T, 2T...20T. 


The program and subroutines are in file STANNARD-P*ANT. 
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APPENDIX D 


APPROXIMATE CALCULATION OF ENHANCEMENT 
OF GAP AREA FOR ELECTROSTATIC ANTENNA 

In construction th« electrostatically - controlled 
membrane minor (ECCM ) , it maybe desiiible to leave a small 
(~lcm) gap at. the point of attachment. If this gap is at 
high potential relative to the nearby election plasma, its 
effective current collection area will be far in excess of 
its nominal area. In this note we provide a rough estimate 
of this enhancement factor. 

A. ESTIMATE OF GAP POTENTIAL 
Consider the geometry of figure 1. 



We will estimate the potential, V^, at the center of the 
gap (width g) by minimizing an approximation to the electro- 
static energy per unit length of gap, U/2 ttR, where 

0 - |E|^ . 

(Minimizing the electrostatic energy is equivalent to solving 
Laplace* s equation) . 

The electric field may be regarded as a superposition 
of the field in the absence of the gap (E= V^/d below the 


m«mbrane; E * 0 above) with the field due to the gap po- 
tential on an otherwise grounded system. To .estimate the 
latter we write the surface potential in the gap region as a 
Fourier sum 

({•Cr) ■ S cos Kjtr. 
it 

Since each component dies exponentially in z: 

<^(r,z) - 2 

the longest wavelength (K^ * ir/g) dominates. Thus we approx- 
imate the gap field by V_/(g/fl*) through a volume 2?rR (g/ir)^. 

The change in electrostatic energy is given by: 



when the first square bracket gives the decrease in energy 
below the membrane and the second the increase above. Mini- 
mizing this expression gives 



Za .JL 

2/ir d 


B. ENHANCED COLLECTION AREA 


The effective collecting surface ^ assumed approximately 
circular) is shown in figure 2. Its height, h, is the distance 



in which the perturbation potential due to the gap is comparable 
to the plasma temperature, 6. Thus 
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-irh/g 


and 


9 - Vq a 


‘aff 


(2irR) (irh). 


Combining thasa rasultS/ wa_hava 

A 


9/ we nave 

2irRg [in (-| jfj)] . 


Note that the enhancement factor (in square brackets) is only 
a logarithmic function of the constants in the problem^ and 
thus not overly sensitive to the approximations used to make 
this estimate. 


Example 


Iii GEO the dominant plasma near a spacecraft is the 
satellite's own photoelectrons, for which 9 leV. If 

V T» SOkeV, g » 1 cm, d » 10 cm, the area enhancement factor 

■ 2 

is 6.7. Since the photocurrent is - 2 x 10 A/ra , the power 

loss for a 100m diameter membrane is: 

Power loss * 6.7 x 2ff x 50 x .01 x (2 x lO"^) x (5 x 10^) 

» 21 watts 
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